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Abstract 

The tunneling through a ferromagnet /normal metal/ferromagnet double junction in the Coulomb 
blockade regime is studied, assuming that the spin relaxation time of electron in the central metallic 
island is sufficiently large. Using the master equation, the current, the tunnel magnetoresistance 
(TMR), and the current noise spectrum have been calculated for devices of different parameters. 
It was shown that the interplay between spin and charge correlations strongly depends on the 
asymmetry of measured device. The charge correlation makes both the chemical potential shift, 
which describes the spin accumulation in the central island, and the TMR oscillated with the same 
period as the Coulomb staircase in current-voltage characteristics. This effect is smeared by the 
temperature. The spin correlation may cause an enhancement of noise at finite frequencies, while 
the zero frequency noise is still always sub-poissonian. The gate voltage causes an oscillation of 
not only conductance, but also TMR and noise. 
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I. INTRODUCTION 



Recently, transport in ferromagnetic single electron transistors (F-SETs) attracts much 
attention from both fundamental and application points of view 111. In such a device a 

n 

small metallic island (or quantum dot (QD) [2j) is coupled via tunneling junctions to two 
leads and capacitively to a gate. Depending on materials of QD and leads (normal (N) or 
'erromagnetic (F) metal) jtoare threejrinds tfF^ET^F/F/F-SET jjflfl, F/F/N-SET 



and F/N/F-SET 



10 



11 



3, Q, 



In these F-SETs the interplay 



between spin and charge correlations during the tunneling process leads to new phenomena 
such as an oscillation of the tunnel magnetoresistance (TMR) with increasing the bias voltage 
[3I or a non-equilibrium spin accumulation on the normal metallic QD [9], which in turn 
may produce a negative differential conductance (NDC) and an accompanied shot noise 



enhancement 



mm 



The aim of this work is systematically to study the spin-dependent transport and the shot 



noise in Coulomb blockade F/N/F-SETs in the sequential tunneling regime [17|, [18]. The 
equivalent circuit diagram of the device is drawn in the inset of Fig. 1(6). We assume that two 
leads are made of the same ferromagnetic metal such as Ni, Fe or Lao^Sro^MnO^, and the 
QD is made of a normal metal with a sufficiently long spin relaxation time such as Al. The 
key parameter in the problem is the polarization defined as P = (Dm — D m )/(DM + D m ), 
where Dm and D m respectively are the density of states for the majority and minority spin 
bands at the Fermi level in leads. Depending on materials, P is often between ~ 0.3(iVi) and 
~ 0.7(La .7Sr 3 MnO s (see [l| p. 149). Such an asymmetry between two spin bands causes 
TMR, defined as TMR = (Rap — Rp)/Rp, where Rp and Rap respectively are device 
resistances in the cases, when magnetizations of two leads are parallel (P) and anti-parallel 
(AP). Neglecting the Coulomb charging effects it was shown that 9j 

TMR = P 2 /(l- P 2 ). (1) 

In order to examine the interplay between spin and charge correlations, by solving the 
master equation, we calculate the current and the TMR for devices different in structure 
parameters, including the polarization P and the lead-magnetization alignment (parallel 
and anti-parallel), at a large range of bias voltage and temperature. The obtained results 
show that the charge effect not only makes TMR oscillated as the bias increases, but also 

n 

considerably enhances it like that observed in F/F/F - and F/F/N-SETs |3j. We also 



calculate the shot noise power as a function of the frequency as well as the bias voltage. 



Though t 



20 



21 



re shot noise has been very extensively studied in metallic QD-structures 



22 



231 ]. for Coulomb blockade F/N/F - single electron tunneling devices, when the 



magnetizations in ferromagnetic leads are collinear, we have found only two works by Bulka 
et al. 
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11] . However, in contrast with our F/N/F-SETs of interest, where, as usually, the 



normal metallic QD is assumed to be large enough so that there is an electron thermalization, 
in these works the QD is assumed to have regularly separated discrete electron levels [h]] or 
to have only a single electron level available for the tunneling process. Such a difference in 
model should be manifested in the noise behavior. Our study shows that though the spin 
correlation may cause a considerable enhancement of noise at finite frequencies, including a 
high peak, the zero-frequency noise is s till always sub-poissonian even in the NDC region. 



The super-poissonian noise observed in [l0|, ITJJ is essentially related to particular structures 
of electron levels modelled. The gate voltage simply leads to an oscillation of the TMR and 
the noise with the same period as that for the conductance. 

The paper is organized as follows. Sec. II is devoted to formulating the problem and 
presenting fundamental expressions. In Sec. Ill we present numerical results of current- 
voltage (I-V) characteristics, of TMR and of noise. In this section the effects of asymmetry, 
of polarization, of temperature, and of gate are in detail discussed. Lastly, a brief summary 
is given in Sec. IV. 



II. GENERAL CONSIDERATION 



Within the framework of the Orthodox theory [18( the state of the F-SET under study 
is entirely determined by the number of excess electrons in QD, n. For a given \n >-state, 
the free energy of the system has the form: 

F(n) = (C 2 V/2-C 1 V/2 + C g V g -ne) 2 /2C t 

- (d7 2 /4 + C 2 V 2 /A + C g V 2 )/2 + e(m - n 2 )V/2, (2) 

where e is the elementary charge, C t = C\ + C 2 + C g , and ni(n 2 ) is the number of electrons 
that have entered the QD from the left (right). Any electron transfer across junctions results 
in a change in free energy F. In the system of interest there are four possible sequential 
electron transfers across two junctions iy = 1 and 2), to the right(+) or the left (-). The 
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change in free energy, Ai 7 ^, associated with these electron transfers depends on the relative 
orientation of magnetizations in ferromagnetic leads. When the two lead-magnetizations are 
in parallel alignment (P-alignment), AF^- is independent of the orientation of electron spin 
[ a —\ (up) or I (down) ] and can be directly defined from (2) as 

AFf (n) = e 2 (1 ± 2n) /2C t T eC g V g /C t T (C 2 + C g /2) eV/C t 
AF 2 ± (n) = e 2 (1 =f 2n) /2C t ± eC g V g /C t =F (d + C fl /2) eV/C t 

When lead-magnetizations are in anti-parallel alignment (AP-alignment), the up (down)- 
spin electron has a larger (smaller) tunneling rate for the source junction than for the drain 
junction. In this case, to ensure a balance between steady tunneling currents for each spin 
orientation, the chemical potential for up (down)-spin state, /^q), in metallic QD should 
be shifted upwards (downwards) by (the same) amount, 8^ = —S/i^ = 5ji. Physically, 
such a chemical potential shift is equivalent to a non-equilibrium spin accumulation in QD. 
Taking into account this chemical potential shift, the change in free energy associated with 
the AP-alignment becomes depending on the orientation of electron spin a and the chemical 
potential shift b\i as : 

AF± (n) = AF±(n) ± (-If" ^ , i/ = 1, 2, (4) 

where AF^(n) is defined in (3) and for the Coulomb blockade device under study, when the 
spin-flip relaxation time is sufficiently large, the chemical potential shift 5/j, can be generally 
calculated from the balanced condition for spin currents. 

Thus, though the free energy F (2) of each charge state is spin independent, due to a 
non-equilibrium spin accumulation in QD and due to a difference between effective tunneling 
resistances associated with the majority and the minority spin, the tunneling rate, in general, 
depends on the spin orientation of tunneling electron. For the electron with spin a, which 
tunnels across the //-junction, to the right (+) or the left (-), the rates are given by: 

T± = (e 2 R ua )^ 1 AF^ T /[exp(AF^ T / k B T) - 1]; a =T, | , (5) 

where AF^ a are corresponding changes in free energy defined in eqs.(3) and (4), R va depends 
on the junction tunneling resistance R u and the spin polarization P as: 

R m) = 2R V /{1 ±P V ), v = 1,2, (6) 

where P\ = P 2 = P for the P-aligment and Pi = —P 2 = P for the AP-alignment device. 
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Using expressions (3)- (6), in principle, we can solve the master equation (ME) or perform 
Monte-Carlo simulations to calculate the current and further the TMR and the noise. In 
practice, for simple structures such as the F-SET under study the ME-method is much more 
efficient. Denoting p(n) as the probability of the state | n > of the system, the ME can be 
written in the matrix form: 

dp(t)/dt = Mp(t), (7) 

where pit) is the column matrix of elements p(n, t) and M is the evolution matrix with 
elements defined as follows: the diagonal elements, 



M(n,n) = -[^r+(n) + ^r- 



ci- 



and non-diagonal ones, 



M(n, m) 



T,a( T t( m ) + r 2 CT ( m )); Hm = n-1 

£ CT (rr CT (m) + r+(m)); if m = n + 1 (8) 

0; otherwise . 



Solving the ME (7) in the condition ^ n p(n, t) = 1, one can calculate the stationary currents 
for each spin channel, I a , o =f or j, which are time- independent and equal to the statistical 
average of corresponding currents through any of junctions, v = 1 or 2: I a =< l vu >= 
e J2nl^ta(. n ) ~ ^va( n )]Pst( n ) > where stationary probability p s t( n ) is defined as p{n <— m\r — > 
°°) = Pst(n)S nm and pin <— m|r) is the conditional probability for having state \n > at 
the time r under the condition that the state was \m > at an earlier time t = 0. This 
conditional probability obeys the same ME as for the probability p(n,t). Then, the total 
stationary net current is simply given by: / = J-f + Jj_. Further, from currents associated 
with the P-alignment and the AP-alignment of lead-magnetizations, we can calculate the 
TMR. 

On the other hand, to study the noise we need to know the time-dependent net current 
/(£), which for F-SETs has the form: 

i(t) = j2ut) = J2J2^(t), (9) 



where I va (t) = e£i[T+ « - r~ t), 9l = C 2 /C t and g 2 = C x /C t [21]- 



In order to calculate the shot noise spectrum S(u) of this current I(t), we extended 



Korotkov's expression of S(u) suggested for metallic SETs 



20| to F/N/F-SETs of interest 



and obtain 

v vuaa' nm 

B n m[^{m\fi')p st {m\fi-) - T^,{m\fi + )p st {m\fi + )]. (10) 

Here, A v = e(I+ + I~) with J± = eJ2 na p st (n)T± a , B = Re[(iul - M) -1 ]; and < > is 

the state obtained from the state \m > by transferring an electron across the //-junction to 
the right (+)/ left (-). 

Thus, once the ME has been solved, one can calculate the stationary current, the TMR 
and the noise. Except the case of zero temperature and at low bias of the first Coulomb 
staircase region, when the current expression can be analytically derived (see Appendix), in 
general, calculations have to be performed numerically. The results obtained for devices of 
different parameters (capacitances, tunneling resistances, polarizations and magnetization 
alignments) in a large range of bias voltage and temperature are presented in the next 
section, where the gate effect is also in detail analyzed. 

III. NUMERICAL RESULTS 

Numerical calculations have been systematically performed for a variety of devices with 
different parameters. The obtained results will be presented for only typical devices with (i) 
symmetric (Ri — R 2 = 16 Rq) and strongly asymmetric tunneling resistances (R± = 1000-Rq 
and i?2 = 5Rq); (ii) two values of polarization, P = 0.3 and 0.9; and (Hi) both P- and AP- 
alignments. Here Rq = h/e 2 is the quantum resistance, and the value P = 0.9 is chosen 
to gain a more clear manifestation of spin effect. In all devices studied the capacitances 
are as follows: C\ = 2C 2 = 2C. With chosen parameters these devices exhibit the most 
interesting behaviors in I-V curves such as a clear staircase structure and an NDC. Regarding 
the elementary charge e, the quantum resistance Rq and the capacitance C as basic units, 
in all figures presented below the voltage, the current, the energy (temperature, chemical 
potential) and the frequency are then measured in units of e/C, e/CRQ, e 2 /C, and (CRq)~ 1 , 
respectively. 

(a) Current and chemical potential shift. Fig.l presents I-V characteristics for two typical 
devices mentioned above, symmetric (a) and asymmetric (6), in a large range of applied 
bias V. In each figure two I-V curves corresponding to the case of P-alignment are certainly 
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coincident (solid lines), independent of P (0.3 or 0.9). In the case of AP-alignment, on 
the contrary, the polarization effect becomes important: the dot-dashed curve of P = 0.9 
is much lower than the correspondingly dashed curve of P = 0.3. In addition, for the 
symmetric device in Fig. 1(a) the difference between two curves, corresponding to the same 
AP-alignment, but with different P (or two curves corresponding to the same P but with 
different alignments) seems to be much larger than that for the asymmetric device in Fig. 1(6) 
and increases with increasing the bias. So, Fig.l shows a relatively stronger spin effect in 
symmetric devices compared to asymmetric ones. This can be understood from the fact that 
the charge effect is more profoundly manifested in the asymmetric device (Fig. 1(6)), where 
I-V curves show clear Coulomb staircases. 




V [e/C] 

FIG. 1: I-V characteristics for (a) the symmetric device (R\ = R 2 = 16i?g) and (6) the asymmetric 
device (i?i = WOORq and R 2 = 5Rq). For both devices: d = 2C 2 = 2C. The currents for P- 
alignment are independent of P (solid lines). For the AP-alignment: dashed line - P = 0.3; 
dot-dashed line - P = 0.9. Insets (top and left in both figures): the chemical potential shift for the 
same device with the same P as in the main figure. The equivalent circuit diagram of the device 
is drawn in the inset in bottom of Fig. 1(6) [without gate and T = 0]. 

For the AP-alignment, when spin relaxation time is sufficiently longer than the time 
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between successive tunneling processes, as mentioned above, a non-equilibrium spin accu- 
mulation in the QD may take place. Such a spin accumulation is equivalent to a difference 
between two chemical potentials, corresponding to two spin states. With increasing the 
bias, /if rises to balance the spin incoming and outgoing rates, and decreases by the same 
amount of Sfi. For systems without Coulomb blockade this chemical potential shift can be 
simply evaluated as Su, = PeV/2 [14], which seems to describe well even the data in the 
inset of Fig. 1(a) for the symmetric device, when the linearity of I-V curves are still weakly 
affected by the charge correlation. For the asymmetric device, as can be seen in the inset of 
Fig. 1(6), though the average value of S/j, does increase with increasing P and V similar to 
that observed in Fig.l(a), the charge correlation makes oscillated with the same period 
of e/2C as the Coulomb staircase in I-V curves. 

In Fig.l the temperature is zero. Fig.2 describes the finite temperature effect in the most 
interesting case of asymmetric devices with P = 0.3 (a) and 0.9 (6). As the temperature 
increases, both the staircase in I-V curves and the oscillation of chemical potential shifts (see 
insets) are steadily smeared. In principle, any finite temperature can destroys the Coulomb 
gap, though in practice the current in the gap may be very small when the temperature is 
still much smaller than the charging energy. Remarkably, for the asymmetric device with 
P = 0.9, shown in Fig.2(6), an NDC region may be appeared at temperatures low enough. 
The observed NDC gradually disappears when the temperature raises. In the particular 
case of appropriate device parameters at zero temperature and low bias, when the current 
expression can be analytically derived as shown in Appendix, actually, we can also derive 
the condition for NDC to be observed. 

(6) Tunnel magnetoresistance. As the typical measure of the spin effect in F-SETs, the 
TMR attracts much attention from both experimental and theoretical sides. For F/N/F- 
SETs under study, as can be seen in Fig. 3, both the value and the bias-dependent behavior 
of TMR are very sensitive to the asymmetry of SETs [(a) and (6) for symmetric, while (c) 
and (d) for asymmetric SETs], to the polarization [(a) and (c) for P = 0.3, while (6) and 
(d) for P = 0.9], and to the temperature [ksT = (solid line), 0.01 (dashed line) and 0.02 
(dot-dashed lines)]. Overall, in agreement with those reported for various F-SETs 3l. la. Il2| 
this figure generally shows that (i) the TMR oscillates with the same period of e/2C as that 
for 5fi and with an oscillation amplitude decreasing as the bias V increases; (ii) the TMR 
generally decreases with increasing V and becomes constant in the large bias limit; and 
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FIG. 2: I-V characteristics for asymmetric devices with P = 0.3 (a) and P = 0.9 (b) at different 
temperatures (ksT): (solid line); 0.01 (dashed line); and 0.02 (dot-dashed line). Insets: the 
chemical potential shift for the same device and at the same temperatures as in the main figure 
[without gate]. 

(Hi) the temperature, smearing the oscillation of TMR at low biases, does not affect TMR 
at large biases, when the Coulomb blockade effect can be neglected. For the symmetric 
F-SETs, (a) and (b), the large-bias limiting values of TMR, defined quite well at eV of 
several charging energy, seem to be in good agreement with the simple expression (1). For 
the asymmetric SETs, (c) and (d), the oscillation of TMR is much stronger and extends to 
a larger range of bias. By comparing four these figures to each other, one can learn much 
about the interplay between spin and charge effects in TMR. In particular, the symmetric 
device with P = 0.9 [Fig.3(6)] shows the most remarkable spin effect with high values of 
TMR, whereas the charge effect is so weak that TMR is almost bias-independent except a 
narrow region close to the Coulomb blockade threshold voltage. 

Takahashi and Maekawa 0] and Kuo and Chen [ljj] have analyzed the gate effect in the 
tunneling current and the TMR in F-SETs. As well-known in the theory of SETs, the gate 
with voltage V g and capacitance C g leads to the effect similar to an offset charge of C g V g , 
which, in turn, can be related to the chemical potential shift in QD. Thus, one can expect to 
see an oscillation of TMR with increasing V g like that observed in the conductance of SETs. 
In Fig. 4 we show the TMR versus the gate parameter C g V g , calculated for the asymmetric 
device with P = 0.3 (a) and 0.9 (b) [the same devices as in Figs. 3(c) and (d)} at two biases: 
V = 0.7 (solid line) and V = 1.4 (dashed lines). It is clear from both figures that all the 
curves, though strongly different in value and form, exhibit a well-defined oscillation with 
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0.20 




V [e/C] 



FIG. 3: TMR is plotted as a function of V [V > V c , where V c is the Coulomb blockade threshold 
voltage] at various temperatures (ksT): (solid line); 0.01 (dashed line); and 0.02 (dot-dashed 
line). Figures: (a) and (b) - symmetric device with P = 0.3 and 0.9, respectively; (c) and (d) - 
asymmetric device with P = 0.3 and 0.9, respectively [without gate]. 

the same period of e/C g as that well-known for the conductance. Note that for the device 
of P = 0.3 (a) the TMR though smaller in value is relatively more sensitive to the change of 
gate, compared to the device of higher polarization in figure (b) . Moreover, comparing two 
curves in each figure, we see that the amplitude of TMR-oscillation is larger at lower bias, 
when the Coulomb blockade effect is more important. 

(c) Shot noise. The shot noise in F-SETs, as mentioned above, has been studied by Bulka 



et al. 
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11], but for QD-models different from ours. Fig. 5 shows the frequency dependence 
of the normalized noise S(u)/2el calculated from eq.(10) for the asymmetric devices with 
AP-alignment. In each figure [(a) for P = 0.3 and (b) for P = 0.9] the results taken at 
different biases, V = 0.47, 0.53, and 0.74, are compared. These biases are specially chosen 
from the typical regions of I-V curve (NDC region, local minimum, and PDC region). It is 
clear in both figures that the normalized noise generally tends to a constant value greater 
than 0.5 in the limit of high frequency. This limiting value depends on device parameters and 
on bias. It may be even greater than 1 as can be seen in Fig. 5 (b). In the opposite region of 
low frequency, the picture is more complicated, very sensitive to device parameters. Even for 
a given device the noise may experience a very high hill or a deep valley, depending only on 
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FIG. 4: TMR is plotted versus the gate parameter C g V g for asymmetric devices with P = 0.3 (a) 
and 0.9 (b) at two biases: V = 0.7 (solid line) and 1.4 (dashes line) [T = 0]. 

the bias. In fact, from the noise expression (10) we can at least show there really exists such 
a maximum (or minimum) of noise at some frequency. The calculation, though elementary, 
is too lengthy to be presented here. We like only to emphasize that the noise enhancement 
observed is due to the spin correlation. This explains why the noise for the device of P = 0.9 
in Fig.5(6) is generally larger than that in Fig. 5(a). As one more demonstration, we show 
in Fig. 5(6) the data for the same asymmetric device of P = 0.9, but with P-alignment, and 
at the same bias V = 0.74 (dashed line). It is clear that in this case like in normal metallic 



SETs 20J there is no noise-enhancement in the whole range of frequency. Here, note that 
for normal metallic SETs the high-frequency limiting value of normalized noise is g\ + of, 
which is equal to 0.5 for the SET with equal junction capacitances \C\ = Cy. 

In Fig. 6 the Fano factor, F n = S(0)/2eI, is plotted as a function of bias V [two upper 
curves, see the left axis] for the same devices as in Fig.5. Interestingly, in both cases under 
study, P = 0.3 (solid line) and P = 0.9 (dashed line), F n oscillates with the same period 
of e/2C as that for S/j,(V) [two lower curves, borrowed from Fig. 2; see the right axis]. An 
accurate comparison of two corresponding curves, F n (upper) and 8[x (lower), for the same 
P reveals that in each period the highest value of F n exactly corresponds to the minimum 
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FIG. 5: The frequency dependence of normalized noise calculated for asymmetric devices with 
AP-alignment and with P = 0.3 (a) and 0.9 (6) at biases (from bottom): V = 0.49, 0.53, and 0.74. 
In (6) the data for the same device of P = 0.9 at V = 0.74, but with P-alignment, is also shown 
by the dashed line, which should be compared with the first solid line in the top [without gate]. 

in 5fi(V). Thus, Fig. 6 demonstrates that the F n - and the 5/i-oscillations have the same root 
as the Coulomb staircase in I-V curves. Experimentally, an exact correspondence between 
the noise oscillation and the current staircase has been observed in the system of a single 



InAs-QD embedded between resonant tunneling barriers [24|. Another important noise 
feature deduced from Fig.6 is that for F/N/F-SETs under study the zero frequency noise 
seems to be always sub-poissonian [F n < 1] in spite of the noise enhancement realized at 
finite frequencies [Fig. 5] and an NDC region observed in I-V curves [Fig. 2]. Deviations of the 
shot noise S(0) from the full (poissonian) value 2eJ, i.e. F n ^ 1, are the subject of a great 
number of works [see references in 23[]. It is widely accepted that the charge correlation may 
suppress or enhance the noise, depending on the conduction regime. And, two phenomena, 
NDC and super-poissonian noise [F n > 1] are not necessarily accompanied by each other 



mm- 



Finally, we demonstrate in Fig. 7 how the gate affects the Fano factor F n for the same 
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FIG. 6: The zero-frequency normalized noises (two upper curves) and the corresponding chemical 
potential shifts (two lower curves) are plotted versus the bias for asymmetric devices with P = 0.3 
(solid line) and 0.9 (dashed line) [without gate]. 

device as in Fig.6(6) [with P = 0.9] and at the same biases [V = 0.49 (solid line), 0.53 
(dashed line), and 0.74 (dot-dashed line)] as in Fig.5. At any bias, as expected, the gate 
makes the factor F n oscillated with the period of e/C g as the gate voltage V g increases. 
The amplitude of oscillation, however, decreases as the bias increases, showing a gradual 
weakness of charge correlation effect. Besides, it is important that, despite the gate-induced 
oscillation, the zero-frequency noise 5(0) is still always sub-poissonian. 

1.0 
0.9 

hH 

Qi 
<N 

X 0.8 
O 

Xil 

0.7 
0.6 

-1.0 -0.5 0.0 0.5 1.0 

C g V g [e] 

FIG. 7: The zero frequency normalized noise is plotted versus the gate parameter C g V g for the 
asymmetric device with P = 0.9 at several biases: V = 0.49 (solid line); 0.53 (dashed line), and 
0.74 (dot-dashed line). 
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IV. CONCLUSION 



We have studied the tunneling through F/N/F-SETs in the Coulomb blockade regime, 
assuming that the spin relaxation time in the normal metallic QD is much larger than other 
characterizing times in the problem. Using the master equation, we have calculated the I-V 
characteristics, the TMR and the current noise spectrum in devices different in structure 
parameters, in polarization, and in lead-magnetization alignment. 

It was shown that the interplay between spin and charge correlations strongly depends 
on the asymmetry of measured device. While in symmetric devices of equal tunneling 
resistances the spin correlation considerably reduces the current, in asymmetric devices 
the charge correlation becomes more important, leading to I-V curves with Coulomb 
staircases and even with NDC regions. The charge correlation also makes both the chemical 
potential shift, which describes the spin accumulation in QD, and the TMR oscillated as 
the bias increases with the same period as the Coulomb staircase. All these effects are 
gradually smeared by increasing the temperature. In the limit of large bias, when the charge 
correlation can be neglected, the TMR becomes independent of bias. Typical manifestations 
of the interplay between two correlations can be also found in the frequency dependence 
as well as the bias-dependence of the noise spectrum. While the charge correlation always 
suppresses the noise in F-SETs, the spin correlation may cause a noise enhancement at 
finite-frequencies, including a very high peak. The zero-frequency noise is however always 
sub-poissonian, regardless of I-V curve behaviors. The gate voltage causes an oscillation of 
not only conductance, but also TMR and noise. While some of current and TMR results 
obtained in this work are generally in agreement with those reported in literature for 
different kinds of F-SETs, the noise results are new and useful to better understand the 
physics of the tunnel process in F/N/F-SETs. 
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Appendix 



Following the way suggested in ref. |25j we can derive an analytical expression for the 
stationary net current I at zero temperature within the first Coulomb staircase region, 
e/2C\ < V < e/Ci (assuming G\ > C 2 and C g = 0). Actually, under these conditions, the 
rate expression (5) becomes simple as Y^ a = 6(— /S.F^ J )\/S.F^ (T \/ e 2 R ua and all the probabili- 
ties p(n) are equal to zero except those for two states | — 1 > and |0 >. The master equation 
can be then exactly solved that gives 

7 _ [r+(-i) + r+(-i)][r 2 + T (o) + r+(o)] 
r+(-i) + r+(-i) + r+(o) + r+(o)- 

Correspondingly, the quantities Ai^(n) in the rate expression are now defined as 

AF+ U) (-1) = -e 2 /2C t - eC 2 V/C t ± // ; AF+ U) (0) = e 2 /2C t - eC x V/C t =F fi , 

where 



li = [eV{l + 7) - y/(l + 7 ) 2 e 2 y 2 - 4(1 - 7 ) 2 «/? ]/2(l - 7), 

7 = (1 - Pf/{1 + P) 2 ; a = edV/Ct - e 2 /2C t ; (3 = eC 2 V/C t + e 2 /2C t ; and C t = G x + C 2 . 
From the current expression obtained the condition for observing an NDC could be also 
derived. 
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